Frequency-dependent current correlation functions from scattering theory 



J. Salo 

Laboratory of Physics, Helsinki University of Technology, PO Box 4100, 02015 TKK, Finland 

F.W.J. Hekking 

\Q . Laboratoire de Physique et Modelisation des Milieux Condenses, 

f^*) ■ CNRS & Universite Joseph Fourier, BP 166, 38042 Grenoble cedex 09, France 

o . 

<N . j.p. Pekola 

Low Temperature Laboratory, Helsinki University of Technology, PO Box 3500, 02015 TKK, Finland 

We present a general formalism based on scattering theory to calculate quantum correlation 
functions involving several time-dependent current operators. A key ingredient is the causality of 
the scattering matrix, which allows one to deal with arbitrary correlation functions. Our formalism 
might be useful in view of recent developments in full counting statistics of charge transfer, where 
detecting schemes have been proposed for measurement of frequency dependent spectra of higher 
moments. Some of these schemes are different from the well-known fictitious spin-detector and 
therefore generally involve calculation of non-Keldysh-contour-ordered correlation functions. As an 
■ illustration of our method we consider various third order correlation functions of current, including 

' the usual third cumulant of current statistics. We investigate the frequency dependence of these 

correlation functions explicitly in the case of energy-independent scattering. The results can easily 
be generalized to the calculation of arbitrary n-th order correlation functions, or to include the effect 
of interactions. 
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PACS numbers: 72.10.Bg, 72.70.+m, 73.23.-b 



I. INTRODUCTION 



Dynamical noise properties of mesoscopic systems have been studied for more than a decade, both theoretically and 
^v^j ' experimentally . By now it is well understood that noise measurements can reveal information on the system that is 
. not contained in its DC conductance. So far, most experiments concentrated on measurement of zero-frequency noise. 
However, several proposals have considered the possibility of detecting finite-frequency noise, for instance through 
emission and absorption measurements using quantum few level systems like quantum dots |2| or small Josephson 
junctions as noise detectors. Successful experiments of this type have been reported recently 0,01 • Finite frequency 
noise is interesting, first of all as one expects the noise to probe the intrinsic dynamics of the conductor and hence the 
noise spectral function should be sensitive to the dwell time td of the carriers. Second, at finite frequency, current 
is no longer spatially homogeneous, and charge piles up in the conductor. Coulomb interaction screens this pile-up 
of charge, at a characteristic charge relaxation frequency 1/r which may well be different from 1/rrj. These issues 
have been studied theoretically for diffusive contacts in Refs. |6|, |jfl . Recent calculations of current noise in chaotic 
cavities 0, Q that take both the energy-dependence of scattering and Coulomb interactions into account show that 
I ■ the frequency-dependent noise spectrum is determined solely by the time r, as long as quantum corrections like weak- 
localization can be ignored. In view of recent interest in the theory of the full counting statistics (FCS) of charge 
transfer ^(ji attention shifted from the conventional noise to the study of the properties of the higher moments. 
Recent measurements have probed the zero- frequency third cumulant [ill IT^ . IT^ | . As far as the frequency dependence 
of the higher cumulants is concerned, the situation changes drastically as compared to conventional noise spectra. 
Calculations of the frequency-dependent third cumulant for a chaotic cavity and for a diffusive conductor 141 show 
marked differences from the conventional noise: it is not only determined by the charge-relaxation time r but also 
jj] ■ shows a low- frequency dispersion that is determined by the dwell time Try. 

A properly designed experiment, capable of measuring the frequency-dependent third cumulant, would thus enable 
one to determine the two relevant time scales separately in a mesoscopic conductor. The question as to how to 
design such an experiment brings us to one of the key problems of this field: what is an adequate detector to 
measure frequency-dependent noise spectra, and which noise spectral function is it actually measuring? Most of 
the applications of FCS discussed so far concentrate on the use of a fictitious spin detector, introduced by Levitov 
and coworkers [TH Il6|. This detector measures Keldysh contour-ordered correlation functions of current. Powerful 
theoretical tools have been developed to calculate these correlation functions; therefore this detector is amenable to 
straightforward analysis. However, the spin detector might not be the most suitable one for detecting finite frequency 
noise. Detectors that interact with the noise source through emission and absorption, like the abovementioned 
quantum detectors might be more suitable for this task. The measured spectra are then not directly related to 
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Keldysh-ordered correlation functions, and different methods are required to determine these spectra theoretically. 

In this paper we develop a method capable of handling arbitrarily ordered correlation functions. The formalism we 
adopt is based on scattering theory [lg , pioneered in 0, 0, ■ It is the natural approach to discuss transport and 
noise in mesoscopic devices. The operator for electric current / is written as the difference between the current carried 
by incident particles /;„ and the current carried by scattered particles I ut : I = I'm — lout ■ The central quantity of the 
scattering approach is the energy-dependent scattering matrix. It must satisfy the causality condition in real-time 
representation, which has immediate consequences for the commutation relations between the operators J; n and J out 
at different times |21[ . As a result, any (anti) time-ordered product of current operators can be conveniently rewritten 
as products of currents I m and I ou t with all in-currents ordered to the right (left) of the out-currents. Denoting 
(anti) time-ordering by T (T) this implies T[/ in (i 1 )i : out (i 2 )] = J out (f 2 )/in(ii) and T[4, (ti)-f ut(i 2 )] = 4(ti)Wt 2 ) 
independent of the ordering of t\ and i 2 ■ This way, the cumbersome time-ordering can be avoided and the remaining 
in-out -ordered products can be readily calculated using the scattering theory. 

We apply the in-out -ordering method to the well-studied case of the third cumulant of charge transfer in a 
mesoscopic conductor. We treat energy-independent scattering, and present the time-dependent cumulant in the 
cases of a tunnel barrier (a quantum point contact), a diffusive wire, and a chaotic cavity. First of all, this enables a 
direct check on the validity of our method. Second, we believe that the zero frequency limit of the calculation provides 
a demonstration of the validity of the result for the third cumulant of a tunnel barrier presented in [22|. This result 
had given rise to some discussion in the literature [Til ITil |23| and methods have been developed to settle the issue in 
a frequency-dependent context |25| . Thirdly, our calculation of the frequency-dependent third cumulant enables 
us to find the asymptotic time-dependence of the third cumulant of the charge transfer, both in the short and the 
long time limits. 

The paper is organized as follows: we first summarize the scattering formalism in order to define the notation used 
later, and use the causality of the scattering matrix to derive important commutation relations between in- and out- 
current operators. They are used to establish operator transformation rules, such as T[Ii n (t±) Jout^)] = Iout(t2)Im(ti), 
which allow one to resolve time-ordered products of currents in terms of in-out -ordered products. Their main 
application is to find multi-current correlation functions, and we explicitly present all three-current correlations, 
which are written in terms of three-current spectral functions of two frequency arguments. To keep the presentation 
transparent, we do not address here issues concerning the finite dwell time of carriers nor do we address interaction 
effects. We thus treat the case of energy-independent scattering where the various spectral functions can be evaluated 
using only the transmission probabilities of the scatterer, valid in the limit where the above-mentioned characteristic 
times td, t vanish. It is important to note that, even though we neglect the energy dependence of the scattering 
matrix, we do respect its causality through the in-out -ordering properties. We finally discuss several different 
detection schemes, which all correspond to different three-current correlation functions and, most importantly, use 
the full-counting statistics approach to derive an expression to the time-dependent third cumulant of transmitted 
charge distribution. 



II. SCATTERING FORMALISM AND CAUSALITY 



A. Scattering theory 

The starting point for our analysis is scattering theory, as developed by Biittiker |l7j| . In this formalism, the current 
operator of non-interacting electrons is given by 

= I E / dEdE'e^-z'^laUmaUE') - bl tn (E)b an (E')]. (1) 

n 

The operators a) an (E) and a an {E) create and annihilate electrons with total energy E in the transverse channel n 
in lead a, incident upon the scatterer. Similarly, the creation b' an (E) and annihilation b an (E) operators refer to 
electrons in the outgoing states. For the two-terminal set-up depicted in Fig. ^ a takes values L and R for the left 
and right leads respectively. The results to be presented below can be easily generalized to any multi-terminal case. 
The creation and annihilation operators obey the anticommutation relations, for instance, 

5 nn ,5(E-E') 

(2) 
0. 

Similar anticommutation relations hold naturally also for operators referring to the outgoing states. 



(E)a an ,(E') + (E)a an (E) = 
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Scatterer 




FIG. 1: Two-terminal scattering problem. Both reservoirs are assumed to be in thermal equilibrium, characterized by a common 
temperature T and potentials such that Vr — Vl — V '. 



The operators a and b are related by the scattering matrix s, 

t>an(E) — 2_j s uf3;nm(E)ap m (E) (3) 

and the creation operators fit and & are correspondingly related by the hermitian conjugated matrix, s^ a a. nm {E) = 

S a/3;mn(E)- 

The matrix s is quite generally unitary and it has dimensions {Nl + Nr) x {Nl + Nr). Its size and the matrix 
elements depend on the total energy E. It has the block structure 



r t' 
t r' 



(4) 



Electron reflection back to the left and right reservoirs is described by the square diagonal blocks r (size Nl x Nl) 
and r' (size Nr x Nr), respectively, while the off-diagonal, rectangular blocks t (size Nr x Nl) and t' (size Nl x Nr) 
determine, in turn, the electron transmission through the sample. 

In order to directly benefit from consequences of causality, we present the current operator as the difference between 
two directed currents, carried by incoming states and outgoing states, respectively [21j. Specializing to the left lead, 
we thus write 



where 



and 



I L {t) = /L,in(*) ~ /L,outW, ( 5 ) 

4in(t) = JE/ dEdE'e^ E - E '^ n al n (E)aLn(E% (6) 



4out(i) = T E / dEdE'e^ E - E '^ h bl n (E)b Ln (E'). (7) 

n 

Now, using Eq. @ as well as its hermitian conjugated version, lL,out(t) can also be written as 

„ /3 _ ^, 7, •/ 



(8) 



a,/3 mnk * 



where indices cv and /3 may take values L or i?. This result makes the dependence of the current operator J out on the 
energy-dependent scattering matrix s(E) explicit. As we will detail below, the commutation properties of directed 
current operators at different times are completely determined by the analytical properties of s(E). 



B. Causality 

In real time, the scattering matrix connects operators of an incoming state with those of an outgoing state by the 
convolution relation 

/oo 
Sap-nm it - T)cL/3 m (r)dr. (9) 
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By causality, the scattering matrix must vanish for negative arguments since otherwise an incident current at t\ would 
cause an outgoing current at t 2 < t\. This is equivalent to requiring that the Fourier transform of the scattering 
matrix, s a p. nm (ui) be analytic in the entire upper half plane, since then 

S a 0;nm{u) = hm / — — -, (10) 

V^0+ J^oo 2m [U)' -U))-IT) 

which can be substituted into the inverse Fourier transform of the scattering matrix in order to obtain 

( \ f did* / f\ f d^ —iujt 1 n{ \ l \ i \ 

S a B;nm\t) = / — — S a f3- nm {L) ) / 7T~ e T~i \ ~ = o(t)S a d- nm {t). (11) 

J-oo 2tt J_ qo 2m [uj 1 -cj)-ir) 

Hence the analytical structure of s(uj) as a function of u> (analyticity in the entire upper half plane) implies causality prl 
127^ . i.e., s(t — t') = if t < t'. Similarly, the hermitian conjugated scattering matrix, ^ a n. nm (uS), must be analytic in 
the entire lower half plane. 

C. Commutation relations 

We will use the analytical structure of the scattering matrix established in the previous subsection, Eq. {UJ, to 
obtain the commutation relations for directed current operators J; n and J ou t at different times |2ll |. Consider the 
commutation relation of Ih,m(ti) and I l ,o\it{i2) ■ Starting from 

[lL,mitl),Ih,oVLt{t2)\ 

= (IT E / dE 1 dE 2 dE 3 dE 4 e i ^- E ^/ h e i ^- E ^/^^ (12) 

III ,«2 

and applying the commutation relations as given in J2J we find that 
[iW*i)A°ut(*2)] 

= (I)" E / dE ld E 2 dE 3 [e^~ E ^e^-^ (13) 

711,712 

_e i ^-^)* I A e i(^-^)t 2 /R^ n2( ^ )sii . ;n2ni(jB3)aini(jE ; 2) ]. 

Integrating over all energies we obtain 

[4,i„(£l),iL,out(i2)] =he 2 ^ 4n 1 (*l) S LL;n 1 n 2 (*2-^l)&L™ 2 (i2)-&L 2 ( i 2)sLL;n 2 n 1 (i2-il)aLn 1 (il)- ( 14 ) 

"l ,« 2 

According to Eq. i|ll|) the commutator l|14f) vanishes identically if t\ is a later instant of time than t 2 [2]]]. We thus 
conclude that 

[h,m(h), i L ,out(t 2 )} cx 6{h ~ h). (15) 

We obtain the commutation relations for lL,m{t\) and -Zl. infe), and for -Z^out^i) and iL,out(^2) using the same 
procedure: both these vanish identically, 

[/L,in(*l),/z,in(*2)]=0 (16) 

and 

[/L,out(*l),4,out(t2)] =0. (17) 

These commutation relations have important consequences for the calculation of time-ordered correlation functions 
involving the operators Ii n (t) and I ut(t), as we will now show. 
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D. Time-ordered correlation functions 

We denote the time-ordering of operators by T[A(t\)Bfa)Cfa) . . .], where the operators appear in descending 
order of times, and the anti-time-ordering by T[A(ti)Bfa)Cfa) . . .], with the opposite order of times. Specifically, 
making use of l|15|). I|l(jfl . and l|17l) . we find the following operator identities: 

T[Lfa)Lfa)} = L(ti)Lfa), 

T[I ou t(tl)I ut(t2)] = lout (tl) lout fa), 

(18) 

T[/in(tl)J ut(t2)] = Ioutfa)Iin(h), 
T[i out fa)Iinfa)] = Iout{h)Iinfa). 

One therefore concludes |2l|: time- ordering a product of directed current operators corresponds to an ordering in which 
all the out-currents I ont are placed to the left of the in-currents I ln . 

As an example, let us consider the two lowest time-ordered correlation functions. Using I(t) — Ii a (t) — I ut(t), one 
obtains 

T[ifa)ifa)] = ii n fa)i in fa) - i utfa)L(h) - iout(h)Lfa) + ioutfa)ioutfa), (i9) 

and 

T[Ifa)ifa)ifa)] = L(tl)iinfa)Lfa) - ioutfa)L(tl)iinfa) - ioutfa)L(h)iinfa) 

- iout(tl)iinfa)iinfa) + ioutfa)ioutfa)iinfa) + iout(h)ioutfa)Lfa) (20) 
+ lout {tl) lout fa) linfa) — Iout(tl)Ioutfa)Iout{t3)- 

For the ordered n-current correlation function, the number of terms containing p out-currents and n — p in-currents 
is just the binomial factor n\/[p\(n — p)\]. The sign of such a term is (— l) p . The anti-time-ordering T can be dealt 
with analogously, but here the in and out currents are ordered oppositely: all the out-currents stand to the right of 
the in-currents. 

The important point here, and one of the central conclusions of Ref. |2l| . is that using in-out ordering one gets rid 
of the cumbersome limits of time integration, normally present in time-ordered expressions. This will enable us in the 
following to straightforwardly calculate Fourier transforms and hence directly obtain the frequency-dependent spectral 
functions of the relevant correlation functions. Moreover, the idea of ordering currents using the in-out formalism is 
quite natural in scattering theory. 

III. IN-OUT THREE-CURRENT SPECTRAL FUNCTIONS 
A. General results 

We now turn to consider various three-current correlation functions of the form (6lL,x{ti)8lL,yfa)o~lL, z fa)), where 
each of the x, y, and z refers to the directed component of the current, either in or out, and SI = I — (I). In the 
time-independent case, they can be expressed using the Fourier transform given by 

(Si L , x fa)Si L , y fa)Si L , z fa)) = J ^V^' 1 "^ J ^ e -^-^S xyz {^ 2 ), (21) 

where S xyz {u>\, u 2 ) are the corresponding three-current spectral functions. (Note that another convention is to take 
the transform with respect to ti — t 2 and t\ — t 3 , which leads to slightly redefined parameterization of the spectral 
functions.) Specializing to the case of equilibrium reservoirs, the spectral functions S xyz (uJi, 0J2) are obtained by 
applying Wick's theorem; we refer the reader to Appendix ^ for details. Specifically, we present results for the 
three-current spectral functions in the general case of an arbitrary energy-dependent scattering matrix in Table [n] 
of Appendix El and for energy-independent scattering in Table IIIII in the same Appendix. Here we just note that 
for the particular case of <Si n ,in,in, the energy integral contains Fermi functions of only one reservoir, and its value 
vanishes then identically This is due to the fact that the in-in-in term does not contain the possibly energy-dependent 
scattering matrix. Spectral functions containing two in-currents also only depend on the Fermi function of the left 
reservoir, but the energy-dependence of the scattering matrix may render the integrals nonzero. Such terms, however, 
vanish in the case of energy-independent scattering so that four spectral functions out of the eight are identically zero. 
The four remaining spectra at zero temperature are depicted in Fig. [3 as functions of the two frequencies u>i and uj 2 . 
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■>! [eV/h] to, iu-| ttj 



FIG. 2: The non-vanishing contributions to zero-temperature three-current spectral functions Si n ,out,out, Sout.in.outi "Sout.out.in, 
and S ut,out,out plotted against the frequencies u>i and U2- The first three of these have values between and ±2Fa = 
±2eV i E„ T n (l - T„) while the last one has the extreme value of -T B = -eV^ J2 n T n (l - T n )(l - 2T„) at the origin and 
it vanishes at high frequencies, |oj2 1 3> eV/h. 



B. Limiting cases of in-out -ordered spectral functions 



Although the true advantage of in-out ordering comes when dealing with general correlation functions, we demon- 
strate here that it also provides a straightforward way to obtain the spectral functions in some special cases which 
have been discussed in literature already earlier. In particular, we investigate here the case of energy-independent 
scattering in the limiting cases in terms of temperature, voltage, and the two frequencies. 

As mentioned above, in the case of energy-independent scattering only four three-current spectral functions out of 
eight possible ones remain nonzero. At zero frequencies, oji — ui2 — 0, only 5 ou t,in,out and <S'out,out,out are finite, with 
their values given by 



S'out.in.out — ~ e V~fo ^""^ T n (l — T n ) — — Ta 

n 

e 3 x - 

<Sout,out,out = — e ^~r / y T„(l — T„)(l — 2T„) = — T± 



(22) 



where Ta = e 2 GVF2 and Tb = e 2 GVFz are expressed in terms of the conductance, G = jJ] n T„, and the Fano 
factors of the second and third order, F 2 = - r„)/£„T„ and F 3 = J2n T n( l ~ T n )(l - 2T„)/^„T„. 

The transmission eigenvalues {T„} are the eigenvalues of the transmission matrix tH. In the high-frequency limit, 
\u>\\, \u>2\ eV/H, the non-vanishing terms are in turn 5i n ,out,out and 5 ut,out,iii) whose values equal 2Fa in the second 
and first octants, respectively, and 5 , ou t l i n ,out = — 2Fa in the first quadrant. 

At finite temperatures such that K\u>2\ *C ksT, the spectral functions become independent of u)\ and oj 2 , and 

the non- vanishing ones are given by 

^in, out, out — 5out,out,in — 

2 3 (23) 

"out,in,out — ^rA- 



IV. DIFFERENT PHYSICAL DETECTOR SCHEMES 



An arbitrary three-current correlation function can always be decomposed into a sum of various in-out -ordered 
spectral functions of the type of Eq. 1)2 whose properties are, at least in principle, known. We will illustrate the 
usefulness of this decomposition scheme now for various examples of three-current correlation functions which have 
appeared in the literature. For simplicity we assume energy-independent scattering such that definite results can be 
obtained for three specific examples. We will first consider accumulated charge by a fictitious spin-detector [l6| , which 
directly depends on the Keldysh-ordered correlation functions, and we use the in-out three-current spectral functions 
to evaluate time-dependent third cumulant of the charge distribution. We also compare this with current statistics 
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derived from an unordered generating function and relate it to some of the results earlier appeared in literature. 
The second example is a classical detector which would correspond to the standard fully symmetrized three-current 
correlation function [2|| , and finally we briefly discuss a partially time-ordered correlation function that appears when 
the time evolution of the density matrix of a multilevel quantum detector is considered, coupled to a non-gaussian 
noise source [2!£ • 



A. Third cumulant of FCS 

The third cumulant of the full-counting statistics, i.e. the first correction term describing the deviation from the 
Gaussian distribution of the charge q transported through the conductor during a time t, has been introduced in 
Refs. |l5j and |l6( , and it is given by 

t i * 

((q 3 )) = I dttf dt 2 J dhS^ih,^,^), (24) 



where the Keldysh-contour ordered correlation function is given by 

s^(h,t 2 ,t 3 ) = \{f[i L {t 1 )i L {t 2 )i L {t z )] + T[i L {t 1 )i L {t 2 )i L {t 5 )] + ^ 

o 

- HhihWLMhfo)) + WL{ti)){i L (h)){i L {h)). 

(25) 

Using the operator relations given by Eqs. (|19[) and H20(l . together with their anti-time-ordered counterparts, and 
regrouping the current operators into deviation operators 5Ii n}OUt (t) = Ii n ,out(t) — (I in, out (t)} , enables one to express 
this particular correlation function as 

S^(t U t 2 ,t 3 ) = (<5/i, in (t 1 )5/L ) i„(t2)5/i ) in(t3) " ^/i 1 in(*l)«ii J in(t2)«/L,out(t3) " ^intM^outfe^info) 

3- - - 3- - 3 - - 

- -^I L ,oxLt{tl)5lL;m{t 2 )8lLAa{h) + ^ ^£,in (O^L,out )<RL,out (^ 3 ) + ^^out (*1 ) Wi.out (fa) in (*3 ) 

- 5lL.out{ti)8I L .out{t 2 )5I L ^t{h))- 

(26) 

Each term here can now be expressed in terms of the Fourier transform of the spectral function, Eq. (|21() . and the 
time integrals of Eq. Ij24(l may be carried out explicitly. This results in 



// 3\\ [ duj i f duj 2 (3) sm(u 2 t) ~ smjuit) + sin[(cji - u 2 )t] 

{{q )} - 2 I ~2tt K (Wl ' W2) oW*!-^) (2?) 



where, for this particular ordering, we have 

3 3 3 

S K (uii,UJ 2 ) = — -5'i n: i njOU t(wi, LU 2 ) — -S'in,out,in( w l7 ^2) ~ 7 *Sout,in,in( w l > ^2) 

3 3 < 28 > 

+ — S'in,out,out('i'lj ^2j + -S'out,out,in(wi, U> 2 ) 

*Sout,out,out 

This result is plotted in Fig. [3] (a) for energy-independent scattering at zero temperature. Note that the multiplier of 
each term in the sum above is obtained with the help of the binomial distribution. The particular ordering for current 
operators, like that in Eq. (|25|l . determines the final weight of each xyz spectral function. 



1. Asymptotic values of the third cumulant 

The third cumulant of FCS can be evaluated in the limits of both short and long times t. For short t the cumulant 
is determined by the values of (uji,ui 2 ) at large frequencies where 5' ou t,out,out(^i, w 2 ) is zero, and 

3 

Sg\u u u2) = 3el/^-^T„(i-T n ) =3r A (29) 

71 
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FIG. 3: (a) Keldysh-ordered, (b) unordered, (c) symmetrized, and (d) non-symmetrized spectra at zero temperature. The 
saturated levels of the spectral functions are all proportional to Fa except for the unordered spectrum, which saturates to zero. 
For the other spectra the level essentially depends on in which areas of the (wi, oj2)-plane the spectrum has non-zero values. 
Both Sin, out, out and Sout.out.m vanish at zero frequency, and the spectral functions (a), (c), and (d) are then determined by 
r B = — Soui, out, out, while the spectrum (b) is given by r B — Ta = 5 ut,in,out — Sout.out.out, as indicated in the graph. 



nearly everywhere in the first quadrant of the {oj\, a->2)-plane and zero elsewhere. Therefore, the short-time value of 
the third cumulant is determined by the Si nj0U t, ut & n d S outj0 ut,in spectral functions since only they have non- vanishing 
high-frequency values. We thus have 

«(,»»« flteV^y;T B (l-r B ) f^i r^ Bin^-siny+sin^-^) 
ft „ Jo 2n J 2n x 1 x 2 {Xi-x 2 ) 

(30) 

= tev^-£r„(i-r n ) = r A t. 

n 

Value of the third cumulant for large t is obtained in a similar manner. As long as (0, 0) ^ 0, the leading order 
is given by 

2tt J 2n u>iUJ2(uJi — UJ2) 



("3") 

For ksT <C eV, only S ou t iO ut,out(0, 0) = — Tb has a non-vanishing value in S K at uj\ = UJ2 = 0, and the linear growth 
at long times is then given by 

«? 3 » = r B t (3i) 

while, in the opposite regime, ksT 3> eV, the directed three-current spectral functions become independent of the 
frequency arguments, and S K = eV\ Y] n T n (l — T n ); the long term cumulant is then given by 

«9 3 )) = r A t. (32) 

Since the Keldysh-ordered spectral function is independent of frequency as long as wi|, 0^2 *C ksT/h, this result 
holds as long as t ^> h/ksT. 

Both these results, Eqs. IKS 11) and l|32ll . are in agreement with those presented in Ref. j^], and thus constitute a 
test of the correctness of our approach. Note in particular that we find ((q 3 ))/t = Tb ~ VL T„,(l — T„)(l — 2T„) fol- 
low temperature. This result has given rise to some discussion in the literature, since Ref. [23 obtained ((q 3 ))/t ~ 
J2 n Tn{l—T n ), different from Eq. I|31(l . Several authors [24l l25| subsequently developed methods to analyze frequency- 
dependent three-current correlation functions in order to assess the correctness of Eq. <|31|) . In Ref. [25j an effective 
action approach together with an involved regularization procedure is used to establish Eq. Ij31|) . According to Ref. [24j 
the frequency dependence of Sk, and hence the result for ((q 3 )), depends on the actual position of the spin-detector 
with respect to the scatterer. Then, both results for ((q 3 )) cited above are found, dependin g on the position of the 
detector. A drawback is that the specific frequency-dependence of Sk postulated in Ref. 2J| generally does not 
conserve current. Let us address the issue here in the framework of the in-out-ordering technique. The Y] — T n ) 

proportionality is obtained in Ref. |23| by considering a straightforward quantum analogue of the classical generating 
function, which leads to the cumulant 



((q 3 )) = (([ drSi(r)^ \ = J dh J dt 2 J dt 3 (sfytflMSIfa)) ■ (33) 
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Tunnel junction, Chaotic cavity Diffusive wire 

0<T<1 P(T) = i^^ p(T) = i 1 

JV \22 TO ^n(l — T n )} T avc g 37; 

^<E„r„(l-T n )(l-2T„)> T avc o M 



TABLE I: Values of the averaged transmission parameters for three different types of noise sources: a tunnel junction, a chaotic 
cavity, and a diffusive wire. Here N is the number of transmission channels, P(T) is the distribution function of transmission 
eigenvalues. In the case of a diffusive wire, L is the length of the wire and I <C L is the mean free path of electrons. 



Note that there is no specific time-ordering in this expression. Use of SI = SI ln — <5J ou t then leads to the entirely 
unordered correlation function 

^Unordered (*1 M, h) = (^iu (* 1 )«L,in (t 2 )«L,m (*3 ) - 6I LM («l )6I L , in (t 2 ) 5I L , out (i 3 ) - 5I LM (t 1 )<5/ L ,out (t 2 )«L,in (t 3 ) 

- <S/L,out(tl)WL,in(£2)WL,in(£3) + ^L,in 1 ) Slh.ont (^) ^L.out (*3 ) + 5I L ,out (t 1 ) <S-/l,,out (*2 ) ^L,in (*3) 
+ WL.out(tl)^L,in(t2)Wi.out(i3) ~ S T L,out (*1 ) Wi.out (*2 ) ^L.out (*3 )) • 

(34) 

The corresponding spectrum is given by 



^Unordered — ~>S'm,in,out i^l , W 2 ) ~ Sin.out.inO^l , W 2 ) — S'out, in.inO^l , W 2 ) /^gs 



~T" ^in, out, out (W1,W 2 ) + S'out ,in,out (W1,W 2 ) + S'out ,out,in (Wl, W 2 ) — S'out ,out,out (wi, w 2 ); 

it is plotted in Fig. |3 (b) for zero temperature. Here two terms on the right hand side of (|35[) contribute at zero 
frequency, namely S ou t,in,out and S ou t,out,out- For the unordered three current correlator, we thus find that the 
corresponding third cumulant is given asymptotically (for large t) by 

{(q 3 )) » (r B - T A )t = -t2eV^ ]Tt„ 2 (1 - T n), (36) 

n 

as found in Ref . j2^| • We therefore conclude that the difference between this result and Eq. 131fl is entirely due to the 
different ordering properties of the two definitions of {(q 3 )), Eqs. (|33|l and l|27|). 



2. Time- dependent third cumulant in various cases 



We consider separately the time-dependent third cumulant generated by three different kinds of noise sources: 
a tunnel junction, a chaotic cavity and a diffusive wire 0, in the limit where intrinsic dynamics and interaction 
effects can be ignored (vanishingly small dwell and charge relaxation times) and scattering can be considered as 
energy-independent. Then, the transmission properties of these noise sources can be summarized as in Table [I] 

In an ideal tunnel junction all the transmission probabilities are small, T n <C 1, and all the three relevant transmis- 
sion quantities are equal, 

T n (l - T„)(l - 2T„) » T »0- -T n )ttJ2 T " = NT °™ ( 37 ) 

n n n 

Here TV is the number of transport modes penetrating the tunnel barrier. Hence, the linear coefficient of the time- 
dependent third cumulant remains the same in both the small and long time limits. Numerical integration of Eq. 
(|27ll demonstrates only this linear increase of the cumulant at all times, as illustrated in Fig. 01(a). 

As can be seen from Table [I] the transmission probabilities of a chaotic cavity on the other hand are symmetrically 
distributed between and 1. Consequently, the coefficient of the out-out-out noise term vanishes and the increase of 
the third cumulant with time is slower than linear, see Fig. 01(b). 

Finally, for a diffusive wire the linear growth dominates again for long times, after an initial transient up to several 
H/eV, as can be seen in Fig.0|(c). 



10 




20 30 

Time [h/eV] 




20 30 
Time [h/eV] 



FIG. 4: The third cumulant for (a) tunnel junction, (b) chaotic cavity and (c) diffusive wire at zero temperature. Both the 
tunnel junction and the diffusive wire show a linear growth at large t due to non- vanishing zero-frequency value of Sout.out.out- 



B. A fully symmetrized three-current correlation function 



A classical noise detector measures essentially a signal proportional to the symmetrized two-current correlation 
function 



1 



(38) 



It is quite plausible to assume that a classical measurement of the third-order correlations would yield a signal propor- 
tional to what is essentially a generalization of QH8JI. i.e., a fully symmetrized three-current correlation function p8l | 



? (3) 
Symin V 



(39) 



- \{m)){m)i(t k )) + i</(ti)></fe)></(t fe )>}. 
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This correlator is indeed symmetric in all permutations of the time arguments t\, t 2 , and £3. We can then immediately 
rewrite the corresponding spectral function with the help of the in-out -ordering technique as 

= 2 '-'in, out, out (^1, id 2 ) + 2*-'iii,°ut,out(w2 — w l) 

+ 2^in,out,out( — w 2, Wi — L0 2 ) + — S ou t,out,ia(~ ^2, ~ ^l) (^) 

+ 2^out,out,in(wi, — U} 2 ) + 2^out,out,in(w2 — k>l,W2) 
— <Sout,out,out (Wl, W 2 ). 

Here the presence of various combinations of u>i and w 2 is due to different orderings of the time arguments t\, t 2 , and 
t 3 , and they also give rise to the hexagonal shape of the spectral function in the (wi, W2) _ plane. This result is plotted 
in Fig. 03(c) which coincides with the one found in Ref. |28j |. 

Comparing Eqs. (|4t)|l and l)28[l. or Figs. El (c) and (a), we see that the symmetrized spectrum is generally quite 
different from the Keldysh contour ordered one. Nevertheless, the two coincide in the zero temperature, zero frequency 
limit such that 5g ymm (0, 0) = Tb and hence corresponds to the usual third cumulant of full counting statistics. 



C. Three-current correlation functions of a multi-level quantum detector 

As it is well-known 0, 0, two- level quantum detector responds to two-current correlators such that the direct 
transition rate to the higher level (absorption), given by the Fermi golden rule, is normally determined by the non- 
symmetrized spectral function 

/oc 
dte iut {SI(t)6I(0)) (41) 
-00 

at the frequency u> — —Aw, where Auj is the level spacing. The corresponding relaxation rate (emission) is given 
by the same spectral function but now at the frequency +Aw. This result can be easily generalized to the case of a 
multilevel detector. 

The next-higher order correction to the transition rate, which includes the effect of transitions via an intermediate 
state of a multi- level detector, depends, among others, on the three-current spectral function Sq '{u)\,u>2), which was 
recently discussed in [2i| l3fj| 

1 /*°° 

S i Q ) (ti,t 2 ,t 3 ) = -— 2 - / duJxd^e-^-^e-^-^S^iuJx,^), (42) 

where the partially time-ordered three-time current correlation function is 

S§\h : t 2 ,t 3 ) = (8I{h)T[SI{i2)5I{ts)]). (43) 

We analyze this correlation function here using the in-out -ordering technique. Expanding in terms of in-out three 
current correlation functions yields 

S§\h,t 2 ,t 3 ) = (S Iin(t 2 )6I ia (t 3 ) - 5I in (h)5I ont (t 3 )5I in (t 2 ) - SI in (h)SI out (t 2 )dI in (t 3 ) 

+ 61^)6 I out (t 2 )SI out (t 3 ) - 6I out {h)5I in {t 2 )6I in {t 3 ) + *Jout(ti)*Wt3)Min(t2) ^ 

+ W«mt(*l)Wout(<2)Win(ta) ~ «out (tl)SI oat (t 2 )SI out (t 3 )) 

such that the corresponding spectral function is 

Sq (u>i,U} 2 ) — <Sin,out,out(^l, ^2) + Sout.out.in (^1 1 — ^2) + 5*out,out,in(wi, LU 2 ) — 5 ou t,out,out (^1, ^2), (45) 

see Fig. |21 (d). The zero temperature, zero frequency limit of this quantity is given by Sq (0, 0) = Tb, i.e., it 
corresponds again to the usual third cumulant of current statistics. 
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D. Discussion 

Apart from the unordered spectral function Eq. 135|) . the various spectral functions discussed so far share many 
common features at zero temperature: (i) None of them contains the Sout.in.out contribution, (ii) The sum of the 

terms containing 0, 1, 2, and 3 out-currents are given by binomial coefficients (— l) k ^ ^ J , where k is the number 

of out-currents. For energy-independent scattering, however, terms with k = 0, 1 vanish, (iii) Regions for which 
|cjJi_2 I > eV/H are only determined by the k = 2 terms (5i ni0 ut,out and S , ou t,out,in) while the zero-frequency value is 
given by the k = 3 term (S'out.out.out)- (iv) In regions where W1.2I > eV/h the value of the spectral function is cither 
zero or it saturates to a constant, unlike the two-current spectrum which increases linearly. The variously ordered 
spectral functions differ mainly from each other based on how the 'spectral power' is distributed in the (wi,^)- 

(3) 

plane: the quantum detector noise Sq has twice the value of the symmetrized noise S , g ymm , but that value is only 
achieved for lo\ > while the symmetrized noise has the constant level everywhere in the (wi,a;2)-plane, except in 
the hexagonal area bound within |wi,2 < eV/h. 

V. CONCLUSIONS 

In this paper we have considered a formalism that facilitates calculation of time-ordered current correlation functions 
and applied it to current noise generated by a phase-coherent scatterer. Causality of the real-time representation of 
the scattering matrix causes products of in- and out-current operators, iin(^i) and loutCfe), to vanish if the in-current 
is taken later than the out-current; consequently, time-ordering of current operators may be expressed using in-out 
ordering, in which the out-current operators stand to the left of the in-currents, and vice versa for anti-time-ordering. 
The in-out ordering can be directly applied to current correlation functions of arbitrary order, and they can be directly 
evaluated in the case of thermal reservoirs. If the scattering matrix is, furthermore, energy-independent the correlation 
functions only depend on the transmission eigenvalues of the scatterer. 

It is highly case-dependent to which particular current correlator a detector responds, and we evaluate three 
alternative functions. While a classical noise detector would respond to a fully symmetrized correlator, the spin 
detector discussed in the case of full counting statistics depends on the Keldysh-contour-ordered correlation function 
and a multi-level noise detector to a partially or fully time-ordered correlator. We obtain all the answers without 
cumbersome time-ordered integrations. 
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APPENDIX A: CALCULATION OF THE THREE-CURRENT SPECTRAL FUNCTIONS WITH 

EQUILIBRIUM RESERVOIRS 

We follow Ref. ^7Jj an d obtain all the three-current spectral functions needed by applying Wick's theorem: 

((4 a * - (4^)) {°Lfln - {a ] p a q - (o+o,))^ 

= (4 a <?)(^ a L)< a n a p) - ( a l a n) (aial)(al r a q ) ( A1 ) 

= SkqSl m S np fk(l — /m)(l — fp) — SknSlpS mq fkfm{^ ~ fp)- 
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Next we insert this result into the expression of a three-current correlation function, such as 
(Si LM {h)SI Liill (t 2 )Si Ltin (t 3 )): 



(6I L , in (ti)SI L , in (t 2 )6I Ltin {t 3 )) 

= |E/ dE 1 dE 2 e i ^- E ^/ h ^J2 J dE * dR 



i(E 3 -E 4 )t 2 /h'_ 



l " " n " p J 

x / (al^E^aLjiEi) - (a^(S 1 )a L)i (^ 2 ))) (al n (E 3 )a L , n (E 4 ) - (aljE 3 )a L , n (E 4 ))^ 
x (al p (E 5 )a L , p (E 6 ) - (al p (E 5 )a L , p (E 6 ))) J 

= (^) 3 ^ J ^ 1 ^3^ 5e i(Bl - £5)tl/ V(^- Bl )* 2 /V^~^)t3/a /i( ^ l)/L(jB3)(1 _ fL{Eb)) 

+ (J)'!] / ^i^3^ 5 e i(El - B3) * l/ ' i e i ( E3 - B5 )* 2 / n e i ( B6 - El ) t3 / ft /i(^i)(l ~ /£(£s))(l - /l(£ 5 )) 
rftJi j[ Ml (ti-t a )+ Bi (t 2 -t 3 }l 



2tt 2tt 



x ^N L JdE f L (E)[l -f L (E + fiwi)][l ~ f L (E + hiu 2 ) - f L (E + hui - hw 2 )\. 



from which we can infer that 



5i„,i„,in(wi, W2) = j dE f L (E)[l ~f L {E + fiw!)][l - + ^ 2 ) - + fo^ - fej 2 )], 



cf. Eq. I|21|l . We make next use of the following results valid for Fermi functions: 

8 Ei 



J dEf(E)[l-f(E + SEi)} = - 
J dEf(E)[l - f(E + SEi)} [1 - f(E + 5E 2 )\ = 1 _ 
J dEf(E)[l - f(E + SEi)]f(E + SE 2 ) = 1 _ ^_ 0SEi 



SE 2 



SE 2 - 6Ei 



I _ e -/3SE 2 i _ e -f3(SE 2 -5E 1 ) 

5E 2 SE 2 - 8Ei 



I _ e -f38E 2 ' i _ e ~p(8E 2 -SE 1 



(A2) 



(A3) 



(A4) 



The integration over energy E can then be performed explicitly. In this particular case of Sin,m,in(wi, u> 2 ), the energy 
integral contains Fermi functions of just one reservoir, and therefore its value vanishes: 



{8lLAn(tl)8I L ; ln (t 2 )5I L , in (t 3 )) = Si at i ni - m (u)i,U 2 ) = 0. 



(A5) 



This is generally true only for <Si n ,min since it does not depend on the possibly energy-dependent scattering ma- 
trix. Spectral functions containing two in currents also have Fermi functions of just the left reservoir, but the 
energy-dependence of the scattering matrix may render the integrals non-zero. Yet in the case of energy-independent 
scattering such spectral functions vanish. 



APPENDIX B: IN-OUT SPECTRAL FUNCTIONS OF THREE CURRENTS 

In Table ITT1 all the eight different three-current spectral functions are listed in the general case of energy-dependent 
scattering and assuming equilibrium reservoirs. The corresponding spectral functions for energy-independent scatter- 
ing are given in Table ITTT1 where {T n } denotes the set of energy-independent eigenvalues of the matrix t't. 



[1] Ya.M. Blanter and M. Biittiker, Phys. Rep. 336, 1 (2000). 

[2] R. Aguado and L.P. Kouwenhoven, Phys. Rev. Lett. 84, 1986 (2000). 



14 



Sin,in,in(wi, <^2) = 
<Jin,in,out (l)1,U) 2 ) 



= -Jd£ ^Tiisl L (E + hu 2 )s L L(E)]f L {E)[l- fLiE + hw!)]^- f L (E + htJ2)] 

-Tr[s{ L (E + huJi)s LL {E + Hcji - huj 2 )]f L {E)[1 - f L (E + hu)i)]f L (E + ho)i - ho> 2 )]\ 

iSin.out^n^l, ^2) 

= - / AE \ti[s\ l {E + Hui)sll(E + Ho) 2 )]f L (E)[l ~ ME + fiwi)][l - f L (E + hu 2 )} 
-Tr[sl L (E + hu>i - huj 2 )s LL (E)]f L (E)[l - f L (E + hui)]f L (E + Hui - hu> 2 )] 

5'out,in,in(Wl,^2) 

= - JdS ^Ft[s\ iL (E)si.i l (E + ?iui)]f L (E)[l - f L (E + licu^l - f L (E + hoj 2 )} 

-Ti[s{ l (E)sll(E + /io;i)]/ L (£)[l - f L (E + hu x )]f h {E + hwi - Huj 2 )]\ 
*J in, out, out (0)1,0)2) 

= (EJ dEiTr[si L (E + fiwi)« ia (B + Hu) 2 )sl a (E + hu: 2 )s LL (E)}f L (E)[l - f L (E + Hu)i)][l - f a (E + ho) 2 )} 

-Tr[s{ L (E + ho>i)s La {E + hu>i - hoi 2 )s\ a {E + hbJi - hcj 2 )s LL (E)]f L (E)[l - f L (E + hui)]f a (E + ho>i - hu> 2 ) 
5*out,in,out (o)i,o> 2 ) 

= lEJdfi {Tr [s La (E + hu)x)sl L {E)s LIj {E + ho) 1 )s{ L (E)] f a (E)(l - f L (E + fiaft))(l - f L (E + ho) 2 )) 
-Tr [s£ £ (-B + hoJi)sl L (E)s Lp (E + nui)sl L {E)^ /l(£)(1 - fg(E + haji))f L (E + Ho)i - hoj 2 )\ 

^outjOutjin 

(o)i,o) 2 ) 

= (EJ d J e|Tr[ s [ L (£) SLa (£ + huj^sl^E + ho)i)s LL {E + Hu 2 )]f L (E)[l - f a (E + tiui)][l -f L (E+ ho) 2 )} 
-Tr[s{ a {E)sLL{E + hui)ai L (E + hu>i - hu) 2 )sL a (E)]f a {E){l - ft (E + hu)i)]f L (E + Hon - hu> 2 ) 

"out, out, out 

= X Zofh I dE { T ^ s L(E)s L ^E + nui)si p (E + h^SL^E + fkJ 2 )s{ 1 (E + hu) 2 )s La {E)] 

xf a (E)[l -f p (E+ hu Jl )][l -f y (E + fiwa)] 
-Tr[sl a (E)s L/3 (E + ho)i)sl (E + hoj 1 )s Lj (E + ho)i - ho) 2 )sl 7 (E + Hon - hu) 2 )s La (E)] 

xf a (E)[l -fp{E + ui)]ME + 0,1 - o) 2 )] J. 
TABLE II: Three-current spectral functions for a general energy-dependent scatterer. 
[3] R.J. Schoelkopf, A. A. Clerk, S.M. Girvin, K.W. Lehnert, and M.H. Devoret, in Ref. pjj. 

[4] R. Deblock, E. Onac, L. Gurevich, and L.P. Kouwenhoven, Science 301, 203 (2003); P.-M. Billangeon, F. Pierre, H. 
Bouchiat, and R. Deblock, Phys. Rev. Lett. 96, 136804 (2006). 

[5] O. Astafiev, Yu.A. Pashkin, Y. Nakamura, T. Yamamoto, and J.S. Tsai, Phys. Rev. Lett. 93, 267007 (2004). 

[6] Y. Naveh, D.V. Averin, and K.K. Likharev, Phys. Rev. Lett. 79, 3482 (1997); Phys. Rev. B 59, 2848 (1999). 

[7] K.E. Nagaev, Phys. Rev. B 57, 4628 (1998). 176804 (2004). 

[8] K.E. Nagaev, S. Pilgram, and M. Buttiker, Phys. Rev. Lett. 92, 176804 (2004). 

[9] F.W.J. Hekking and J. P. Pekola, Phys. Rev. Lett. 96, 056603 (2006). 
[10] Quantum Noise in Mesoscopic Physics, edited by Yu. V. Nazarov (Kluwer, Dordrecht, 2003). 
[11] B. Reulet, J. Senzier and D. Prober, Phys. Rev. Lett. 91, 196601 (2003). 

[12] Yu. Bomze, G. Gershon, D. Shovkun, L.S. Levitov, and M. Reznikov, Phys. Rev. Lett. 95, 176601 (2005). 

[13] S. Gustavsson, R. Leturcq, B. Simovic, R. Schleser, T. Ihn, P. Studerus, K. Ensslin, D. C. Driscoll, and A. C. Gossard, 

Phys. Rev. Lett. 96, 076605 (2006). 
[14] S. Pilgram, K. E. Nagaev, and M. Buttiker, Phys. Rev. B 70, 045304 (2004). 
[15] L.S. Levitov and G.B. Lesovik, cond-mat/9401004 (unpublished). 
[16] L.S. Levitov, H. Lee, and G.B. Lesovik, J. Math. Phys. 37, 4845 (1996). 
[17] M. Buttiker, Phys. Rev. B 46, 12485 (1992). 

[18] G.B. Lesovik, Pis'ma Zh. Eksp. Teor. Fiz. 49, 513 (1989) [JETP Lett. 49, 592 (1989)]. 
[19] B. Yurke and G.P. Kochanski, Phys. Rev. B 41, 8184 (1990). 



15 



5j»,i»,i»(wi)k>2) =0 
Si»,in,out(wi,a;2) =0 
Shi, out, in (wi, W2j =0 
Sout,in,ln(wi,o;2) =0 

5 ln ,out,out(wi, wa) = £ E„ T n (l - T n ) J dEU L (E)[l - f L (E + tiui)][l - f R {E + hu> 2 )} 

-f L (E)[l - f L (E + hu3!)]f R (E + twi - feja)]| 

5out, to ,out(wi, w 2 ) T "( x - T «) I dEU R (E)[l - Sl(E + Hui)][l -f L (E + fkj 2 )} 

-!l{E)[1 - f R (E + Hui)]f L (E + hui - hto 2 ) 

S'out,out,m(a;i,W2) =xEn T «( 1 - T ») / dEU L {E)[l - f R (E + hui)][l - f L (E + hcj a )] 

-f R (E)[l - f L (E + HuJi)]f L (E + hui - Huj 2 )] 

S , »t, ut,o U t(wi,Wa) = £ E„ r„(l - T n ) 2 / dB|/i(B)[l - f L {E + hui)][l - f R (E + hu 2 ) - f R {E + fax - huj 2 )] 

+f L (E)[l - f R (E + hui)][l - f L {E + hw 2 ) - f L (E + Hun - hw 2 )} 
+f R (E)[l - f L (E + hu!)][l -f L (E + fkj 2 ) - f L (E + hui ~ huj 2 )\ 



+ 17 En T n(! - T ") / dEU L {E)[l - f R (E + fiwi)][l - + Bw 2 ) - M-E + - fe^)] 

+f R (E)[l -f L (E + hui)][l - f R {E + huj 2 ) - f R (E + fvui - hto 2 )] 
+f R (E)[l - f R (E + hu>i)][l — fi, (E + hw 2 ) - f L (E + hui - huj 2 )] 

TABLE III: Three-current spectral functions for an energy-independent scatterer. 



[20] M. Biittiker, Phys. Rev. Lett. 65, 2901 (1990). 

[21] C.W.J. Beenakker and H. Schomerus, Phys. Rev. Lett. 86, 700 (2001). 
[22] L.S. Levitov and M. Reznikov, Phys. Rev. B 70, 115305 (2004). 

[23] L.S. Levitov and G.B. Lesovik, Pis'ma Zh. Eksp. Teor. Fiz. 55, 534 (1992) [JETP Lett. 55, 555 (1992)]. 
[24] G.B. Lesovik and N.M. Chtchelkatchev, Pis'ma Zh. Eksp. Teor. Fiz. 77, 464 (2003) [JETP Lett. 77, 393 (2003)]. 
[25] A.V. Galaktionov, D.S. Golubev, and A.D. Zaikin, Phys. Rev. B 68, 085317 (2003); Phys. Rev. B 68, 235333 (2003). 
[26] K. Shepard, Phys. Rev. B 43, 11623 (1991). 

[27] M. Biittiker, A. Pretre, and H. Thomas, Phys. Rev. Lett. 70, 4114 (1993). 

[28] D. Golubev, A.V. Galaktionov, and A.D. Zaikin, Phys. Rev. B 72, 205417 (2005). 

[29] T. Ojanen and T. T. Heikkila, Phys. Rev. B 73, 020501(R) (2006). 

[30] V. Brosco, R. Fazio, F.W.J. Hekking, and J.P Pekola, cond-mat/0603844 (unpublished). 



